Non-Markovian reduced dynamics and entanglement evolution of two coupled spins in 

a quantum spin environment 
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The exact quantum dynamics of the reduced density matrix of two coupled spin qubits in a 
quantum Heisenberg XY spin star environment in the thermodynamic limit at arbitrarily finite 
temperatures is obtained using a novel operator technique. In this approach, the transformed 
Hamiltonian becomes effectively Jaynes-Cumming like and thus the analysis is also relevant to cavity 
quantum electrodynamics. This special operator technique is mathematically simple and physically 
clear, and allows us to treat systems and environments that could all be strongly coupled mutually 
and internally. To study their entanglement evolution, the concurrence of the reduced density 
matrix of the two coupled central spins is also obtained exactly. It is shown that the dynamics 
of the entanglement depends on the initial state of the system and the coupling strength between 
the two coupled central spins, the thermal temperature of the spin environment and the interaction 
between the constituents of the spin environment. We also investigate the effect of detuning which 
in our model can be controlled by the strength of a locally applied external magnetic field. It is 
found that the detuning has a significant effect on the entanglement generation between the two 
spin qubits. 

PACS numbers: 75.10.Jm, 03.65.Yz, 03.67.-a 



I. INTRODUCTION 

One of the most promising candidates for quantum 
computation is spin systemsii 2 -^^^^ due to their 
long decoherence and relaxation time . Combined with 
nanostructure technology, they have the potential ad- 
vantage to scale up to large systems. Just as other 
quantum systems, the spin systems are inevitably influ- 
enced by their environment, especially the spin environ- 
ment. As a result, decoherence due to the presence of 
the environment will cause the transition of a system 
from pure quantum states to mixed ones. The decoher- 
ent behavior of a single spin or several spins interacting 
with a spin bath has attracted much attention in recent 
year o 10 ' 11 ' 12 ' 13 . The interaction in such a spin bath sys- 
tem often leads to strong non-Markovian behavior. The 
usual Markovian quantum master equations which are 
widely used in the area of atomic physics and quantum 
optics may fail for many spin bath models. Therefore 
it becomes more and more important to develop meth- 
ods that are capable of going beyond the Markovian 
approximatio n 14 ' 15 ! 16 . 

Entanglement has been recognized as one of the most 
amazing aspects of quantum mechanics. It has been con- 
sidered as important resources for applications in quan- 
tum communication and information processing, such 
as quantum teleportation 17 , quantum cryptography 18 , 
quantum dense coding 1 ^, and telecloning 20 . It is also 
believed to be one of the features that make quantum 
computers more powerful than classical ones. For spin 
systems, much attention has been dedicated to the prob- 
lem of thermal entanglemen t 21 : 22 ' 23 , i.e., to quantify en- 
tanglement arising in spin chains at thermal equilibrium 



with an environment or a reservoir. In this approach, a 
thermal distribution of the system energy levels is deter- 
mined by the environment temperature, but the detailed 
interaction between the system and environment, and the 
evolution of the system toward the thermal equilibrium 
are explicitly ignored. 

A quantum system exposed to environmental modes is 
described by the reduced density matrix when the envi- 
ronment modes are traced over. The time evolution of the 
reduced density matrix is usually very difficulty to obtain 
in the case of non-Markovian process. Recently, the dy- 
namics of the reduced density matrix for one-, two- and 
three-spin-qubit systems in a spin bath described by the 
transverse Ising model has been analyzed without making 
the Markovian approximation, but using a perturbative 
expansion method 2 ^ or a mean-field approximatio n 25 : 26 . 
The interaction between the system and the spin bath for 
these cases was assumed to be of a Ising type. It has also 
been reported recently that the exact reduced dynam- 
ics and for one- and two-spin-qubit systems in a spin- 
star environment^ 7 , has been derived and analyze d 14 ' 15 . 
There, the interaction between the system and environ- 
ment was assumed to be of a Heisenberg XY interaction. 

In this paper, we study a two-spin-qubit system in a 
spin star configuration, similar to the case studied in 
Ref. [l6|. There are however several important differences 
between our model and that of Ref. [l6|. In Ref. [lj| the 
interaction between the two-spin qubits and the inter- 
action between the constituents of the spin environment 
are neglected, i.e., no internal dynamics for both the spin 
qubit system and the spin environment is considered. In 
addition, the spin environment is assumed to be initially 
in an unpolarized infinite temperature state. As a re- 
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suit, no dependence of the environment temperature on 
the dynamics and entanglement is present. It is under 
these conditions that the "exact" dynamics is reported 
in Ref. 16. Neglecting the direct interaction among the 
constituents in the environment and considering only the 
infinite temperature initial state may not be proper in 
dealing with spin baths. In this paper, we investigate 
a more general case. Using a novel operator technique, 
we present an exact calculation of the dynamics of the 
reduced density matrix of two coupled spins interacting 
with a thermal spin bath at finite temperatures in the 
thermodynamic limit. In our model, the interaction be- 
tween the constituents of the spin environment, the inter- 
action between the two spin qubits, and the interaction 
between the the spin qubit system and the spin environ- 
ment are all of the Heisenberg XY type and can all be 
taken into account simultaneously. In addition, we in- 
clude also the Zeeman coupling between the spin qubits 
and a locally applied external magnetic field. To quan- 
tify quantum entanglement dynamics of the two coupled 
spins under the influence of the spin bath at arbitrar- 
ily finite temperatures in the thermodynamics limit, we 
calculate the exact time evolution of their concurrence. 
Our model involves the Heisenberg XY coupling which 
has extensive applications for various quantum informa- 
tion processing proposal o 28 ' 29 ' 30 ' 31 ' 32 ' 33 ' 34 ' 35 . In addition, 
the transformed Hamiltonian of the total system in our 
approach becomes effectively Jaynes-Cumming like and 
thus our analysis is also very relevant to cavity quantum 
electrodynamic s] 28 ' 29 ' 34 ' 35 . 

The paper is organized as follows. In Sec. II the model 
Hamiltonian is introduced and the operator technique is 
employed to obtain the reduced density matrix, taking 
into account the memory effect of the environment. From 
the reduced density matrix, the entanglement measure of 
concurrence of the coupled spin system is calculated in 
Sec. III. Conclusions are given in Sec. IV. 



II. MODEL AND CALCULATIONS 

We consider a two-spin-qubit system interacting with 
bath spins via a Heisenberg XY interaction. The system 
and bath are composed of spin- 1 ; atoms. We restrict our- 
selves to a star-like configuration with couplin g of equal 
strength, similar to the cases considered in Refs. [l4]il5li27l . 
The interactions between bath spins are also of XY type. 
In Refs. [2lf25l[26l . a similar but somewhat different type 
of Ising interactions between the constituents of the spin 
bath was considered. The Hamiltonian for the total sys- 
tem is 



H = He 



H 



SB 



Hp 



(1) 



Here, H$ and Hb are the Hamiltonians of the system and 
bath respectively, and H$b is the interaction between 
the m 15 ' 27 . They can be written as 



He 



MoOSqi + Sq 2 ) + (S \S 02 + S 01 Sq 2 ) 



(2) 



H 



SB 
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N 



N 



i=l 



(4) 



i¥=3 



where fj,o represents the coupling constant between a lo- 
cally applied external magnetic field in the z direction 
and the spin qubit system. f2 is the the coupling con- 
stant between two qubit spins. S^ and S^ (i=l,2) are 
the spin-flip operators of the qubit system spins, respec- 
tively, and S~ are the corresponding operators of 
the ith atom spin in the bath. The indices of the sums 
for the spin bath run from 1 to N, where N is the num- 
ber of the bath atoms, go is the coupling constant be- 
tween the qubit system spins and bath spins, whereas g 
is that between the bath spins. Both coupling strengths 
are rescaled such that the free energy is extensive and a 
nontrivial finite limit of N — > 00 exist s 15 ' 36 . 

By using collective angular momentum operators J± = 
J2i=i we rewrite the Hamiltonians, Eqs. (3) and (4), 
as 



H SB = 

H B = 



9o 



[(#1 



S02) J- 



(^01 



2j 



(J+J-+J-J^ 



9, 



J+] . (5) 
(6) 



where j — N/2 is the length of the pseudo-spin. After 
the Holstein-Primakoff transformation 3 - 7 -, 

J+ = b+ (y/2j - b+b) , J_ = [y/2j - b+b) b, (7) 

with [b,b + ] = 1, the Hamiltonian, Eqs. (5) and (6), can 
be written as 



H 



sb = go 



H B = g 



(S + i + S+) \p^-b + (S 1 + S 1) b + \jl 



N 



TV V N 



In the thermodynamic limit (i.e. 
temperatures, we then have 



H SB 
H b 



9o [{So\ 
2gb + b. 



N 



[Sqi + S) 



00) at finite 

1. (10) 

(11) 



Equations (2), (10) and (11) are then effectively equiva- 
lent to the Hamiltonian of a Jaynes-Cumming type. They 
describe two coupled qubits interacting with a single- 
mode thermal bosonic bath field, so the analysis of the 
problem is also relevant to cavity quantum electrodynam- 
ics quantum information processing proposal s 28 ' 29 ' 34 ' 35 . 
We note here that due to the high symmetry of our 
model, the coupling to the environment is actually rep- 
resented by a coupling to a single collective environment 
spin. After the Holstein-Primakoff transformation and 



b+b 



-9- (9) 
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in the thermodynamic limit, this collective environment 
spin is transformed into a single-mode bosonic thermal 
field. The effect of this single-mode environment on the 
dynamics of the two coupled qubits is extremely non- 
Markovian. This reflects onto, for example, the revival 
behavior of the reduced density matrix or entanglement 
evolution of the two coupled spins, which will be shown 
later. This is different from the usual environment mod- 
els which consist of very large degrees of freedom (e.g. 
many bosonic modes) and often cause the reduced dy- 
namics of the system of interest displaying an exponential 
decay in time behavior. So the Markovian approxima- 
tion usually used in quantum optics master equation will 
not work in our model. One may perform perturbation 
theory for weak-coupling case, but the single-mode envi- 
ronment in our model will not remain in thermal equilib- 
rium state as usually assumed for an environment with 
very large degrees of freedom in the weak-coupling mas- 
ter equation approach. 

Using a special operator technique, we can obtain the 
exact reduced density matrix for the two coupled qubits 
by tracing over the degrees of freedom of the bosonic bath 
at arbitrarily finite temperatures. Reference l38l reported 
the theoretical results of entanglement dynamics of a cou- 
pled two-level atoms interacting with a cavity mode em- 
bedded in an effective atomic environment. However the 
influence of the environmental temperature was not con- 
sidered. In Ref. 39, a decoupled two-qubit system inter- 
acting with a single-mode thermal field at resonance (i.e., 
zero detuning) in the context of cavity electrodynamics 
was studied. There, the dynamics of the reduced density 
matrix for the two-qubit system is obtained using the 
method of the Kraus operator representation. In this pa- 
per, we use a different approach of operator technique to 
obtain the exact non-Markovian dynamics of the reduced 
density matrix for the two-qubit system for our model of 
Eqs. (2)-(4) with the bath spins in the thermodynamics 
limit, or equivalently Eqs. (2), (10) and (11). Differ- 
ent from the case considered in Ref. 39, our model fur- 
thermore includes the coupling between the two qubits 
and investigates the effect of detuning (i.e., the single 
bosonic bath mode is not necessarily resonant with the 
qubit transition frequency). The detuning in our model 
is represented by (/xq — 2g) and it could be controlled by 
the strength of a locally applied magnetic field, i.e., the 
po term in Eq. (2). We find that the detuning has a sig- 
nificant effect on the entanglement generation between 
the two qubits. 

We assume the initial density matrix of the total sys- 
tem to be separable, i.e., p(0) = \ip)(ip\®PB- The density 
matrix of the spin bath satisfies the Boltzmann distribu- 
tion, that is ps = e~ HB / T /Z, where Z = Tr (e~ HB / T ) is 
the partition function, and the Boltzmann constant has 
been set to one. At absolute zero temperature, no exci- 
tation will exist. The bath is in a thoroughly polarized 
state with all spins down. With the increase of tempera- 
tures, the number of spin up atoms increases. Note that 
in Ref. HtH the non-interacting bath spins are assumed to 



be initially in the unpolarized infinite temperature state. 
The most general form of an initial pure state of the two- 
qubit system is 



|^(0)) = a|00) + /3|11) + 7|01) + <J|10>, 
with - 1 1 



i. 



(12) 
(13) 



We might proceed the calculation with this general initial 
state, but the final analytical solution would, however, be 
somewhat complicated. For analytical simplicity, in the 
following we set 7 = 6 = 0. We note that the general 
initial qubit state case can be calculated in a similar way 
presented below. 

By taking the initial state of the two qubit system to 
be |^) = a 1 00) + /3|11), the reduced density matrix can 
be written as 



Ps(t) = |M 2 tr B 
+-ap*tx B 
+ la*{3tr B 



-iHt 



|00>€ 



'<00|e 



-iHt 



-iHt 



|00)e- ffB / T (ll|e <ff * 



ll)e- ffB / T (00|e' 



where 



--1/91^ 



z = 



,-iHt ln)e -H B /T (nle i 



1 



Hi 



HI 



,-2g/T' 



(14) 



(15) 



The matrix p s (t) is a 4 x 4 matrix in the standard basis 
1 00), 1 1 ) , 1 10) , 1 1 1) . In order to obtain the exact reduced 
density matrix elements, we have to evaluate Eq. (14) ex- 
actly. However, it is difficult to do so using usual meth- 
ods, because the treated system and environment could 
all be strongly coupled mutually and internally. In the 
following, we will present a special operator technique 
to obtain the exact density matrix elements. As shown 
below, our treatment is mathematically simple and phys- 
ically clear, and may be easily extended to more compli- 
cated systems with strong coupling. Note also that our 
method also applies to the case that the two-qubit system 
is initially in a mixed state. For example, if the initial 
state for the qubits is p s (0) = |a| 2 |00)(00| + |/?| 2 |11)(11|, 
the corresponding reduced density matrix is Eq. (14) pro- 
vided that the second and third terms on its right hand 
side are removed. 

The basic idea of our operator technique is as follows. 
Before tracing over the environmental degrees of freedom, 
we will first convert the time evolution equation of the 
qubit system under the action of the total Hamiltonian 
into a set of coupled non-commuting operator variable 
equations. Then by introducing a new set of transforma- 
tion on the operator variables, we turn the coupled non- 
commuting operator variable equations into commuting 
ones. As a result, they can be solved exactly by using 
the general method of solving coupled first-order ordi- 
nary differential equations for ordinary variables. After 
that, the trace over the environmental degrees of freedom 
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can be performed and the exact reduced dynamics of the 
qubit system can be obtained. 

From the total Hamiltonian H, we can see that it con- 
sists of operators b, b + , and S^, where S^, and 
change the ith (i = 1, 2) qubit spin from state |l)i to |0)i, 
and vice versa. It is then obvious that we can write in a 
most general form that 



-iHt 



|11) = A\00) + B\01) + C\W) + D\ll), (16) 



where A, B, C, and D are functions of operators b, b + , 
and time t. Using the Schrodinger equation identity 



• « (g-iHtjix)) = H{e- iHt \ll)) 



dt 

and Eq. (16), we obtain 



(17) 



u = 

dt 


-i [-fi A 


4- 2gb + bA 4 


-gob+B + gob+C] (18) 


dt 


-i [g bA + 


2gb+bB + 


QC + g Q b+D] , (19) 


dt 


-* [9ob A + 


ilB + 2gb 


f bC + g Q b + D] , (20) 


d 

—D = 

dt 


-i [ 9o bB -f 


- gobC + no 


D + 2gb + bD] , (21) 



with initial conditions from Eq. (16) being A(0) = 0, 
5(0) = 0, C(0) = 0, and D{0) = 1. As A, B, C, D 
are functions of b + and b, they are operators and do not 
commute with each other. Equations (18)-(21) are thus 
coupled differential equations of non-commuting operator 
variables, which can not be solved by using conventional 
methods for ordinary number variables. 

The crucial observation to solve the problem is that 
the Hamiltonian, Eqs. (2), (10) and (11), is of an effective 
Jaynes-Cumming type and it can be block-diagonalizcd 
in the dressed state subspace of |i, j; n), with i + j + n = 
constant. Here \i,j) represent the qubit states and \n) 
are the bosonic field number states. As a result, we may 
rewrite Eqs. (18)-(21) in such a subspace. By introducing 
the following transformation 



A = b+b+e-^(b + b+i)t Au 

B = b+er a ^ b+b +^B 1: 

C = b+e-^f+^d, 

D = e -iMb + W)t Dl: 



equations (18)-(21) then become 
d 



dt 
d 

dt 
d 

dt 
d 

Tit 



At 



B t 



Ci 



i(/x - 2g)At - ig {Bt + Ct) , 
-i [.go (2 + h)At + flCt + goD{\ 
-i [.go (2 + n)At + flBt + g D{\ 



(22) 
(23) 
(24) 
(25) 



(26) 
(27) 
(28) 



where n = b + b. Note that for initial qubit state in |11) 
on the left hand side of Eq. (16), the transformation, 
Eqs. (22)-(25), is chosen in such a way that the bosonic 
field operator(s) in front of the exponential term in 
Eqs. (22)- (25) together with its corresponding qubit state 
on the right hand side of Eq. (16) make i + j + n a con- 
stant value and the initial condition D(0) = -Di(O) = 1. 
That is, the operator coefficient A of |00) state on the 
right hand side of Eq. (16) requires two field creation 
operators in Eq. (22), B and C of 1 01) and 1 10} respec- 
tively require only one field creation operator in Eqs. (23) 
and (24), and D of |11) does not need any field opera- 
tor in Eq. (25). The exponential term in Eqs. (22)-(25) 
is introduced to make the resultant equations more con- 
cise. As a consequence, the coefficients of Eqs. (26)-(29) 
after the transformation (22)- (25) involve only the oper- 
ator n. Therefore At, Bt, Ct, and Dt are functions of 
n and t, and commute with each other. We can then 
treat Eqs. (26)-(29) as coupled complex-number differen- 
tial equations and solve them in a usual way. This novel 
operators approach thus allows us to solve Eq. (16) and 
then consequently the non-Markovian dynamics of the 
reduced density matrix of the qubit system. 

We note again that the crucial point of the method 
used here is to find proper transformations to change the 
coupled differential equations of non-commuting operator 
variables to the coupled differential equations of complex- 
number variables. This can be done when the effective 
Hamiltonian can be block-diagonalized. This is the case 
of Jaynes-Cumming model and other models which con- 
tain interaction Hamiltonian of the forms of, for example, 
(S+ + S+) bb + (S^t + S02) b + b+, S+S+b + SoiS - 2 6+, 
etc., regardless how strong their interaction strengths 
are. If the total effective Hamiltonian can not be block- 
diagonalized, for example, for the effective spin-boson 
model in Ref. 0, the operator method used here will 
then not apply to solve the problem exactly. 

As we are working in the Schrodinger picture, the basic 
operator n = b + b is time independent (sometimes the 
operators could have time dependence explicitly; however 
this is not the case here). From Eq. (16) and Eqs. (22)- 
(25), the initial conditions at t = are given by 



At(Q) 
Bt(0) 
Ct(0) 
Dt(0) 



(30) 

(31) 
(32) 
(33) 



-ig (l + ft) (Bt + Ct) - i(fx ~ 2g)Dt,(29) 



In general, we can solve Eqs. (26)-(29) exactly via the 
initial conditions Eqs. (30)- (33). As we aim to obtain 
analytical expressions for the reduced qubit dynamics, 
for the sake of analytical simplicity we consider the on- 
resonant case, i.e., ^0 = 2g. We can easily tune the lo- 
cally applied external magnetic field to satisfy this con- 
dition. We will give the numerical results for the off- 
resonant case in Fig. 6. We then obtain for the on- 
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resonant case 

M - 1 



+ 



Bi = C x = - 



3 + 2ft y/fli + 8(3 + 2n) 3 2 
5o 



Ai A2 



04) 



v/ft 2 + 8(3 + 2n)<? 2 
2 + h 2 ff 2 (l + n) 



"giAit giA 2 t 



3 + 2n ^! 2 +8(3 + 2ft) 5 2 



Ai A 



, (35) 



where 



Al,2 



-fi± ^/fi 2 + 8(3 + 2n) ff 2 



(37) 



Following the similar calculations above, we can evaluate 
the time evolution for the initial two-qubit spin state of 
1 00). Let 

e - im |00) = £?|00) + F|01) + G|10) + Kill). (38) 
In a similar way, we have 

E = e- i2g ( b+b - l >E 1 , 



F = be-M^-^Fi, 
G = be-M^^d, 
K = bbe-M b+b -^ Kl , 



(39) 
(40) 
(41) 
(42) 



and then obtain 
n — 1 

Ei = 



2gln 



F\ = G\ = — 

K x = 



.9o 



Ai A 



y/Q? + 8 (2n - l)gl 
-1 2 5o 2 



giAit _ giA 2 t 



(44) 



2n - 1 v /r! 2 + 8(2n-l) ffo 2 



giA'jt giAit 



where 



-n± + 8(2n- 1) 5 2 



(46) 



From Eq. (14) and all the results that we obtained, the 
reduced density matrix can be written as 



Ps(t) = 



/ pn P14 

/022 P23 

/O32 P33 

\ p* 4 ,044 



(47) 



where 



OO 



+ e -2 S r l /T 



n=0 



f± £ AiA+(n + 1)(„ + 2)e- 2 ^ T ,(48) 



Pl4 = afr^^D+e-We"*, 

P22 = P23 = P32 = P33 

00 

= \ a \^J2 F ^ ne ~ 29n/T 
n=l 

+ |/3| 2 i£s 1 B+(n + l)e- 2 WT ; 



(49) 



(50) 



n=0 
00 



P44 



- ]T KiKfn(n - l)e 

n=2 
-. 00 



2gn/T 



,-2gn/T 



(51) 



n=0 



In Eqs. (48)-(51), the trace over the environmental de- 
grees of freedom has been performed and the operator 
n has been replaced by its eigenvalue n. In a similar 
way, the solutions for the reduced dynamics of the two- 
coupled spins under the influence of the quantum Hciscn- 
berg XY spin star bath in the thermal dynamics limit at 
arbitrarily finite temperatures for arbitrary initial states 
of |V(0)> = a|00)+/3|H)+7|01)+<J|10) can be obtained. 



III. CONCURRENCE AND ENTANGLEMENT 
DYNAMICS 

We use the concurrence^ 1 - to measure the entanglement 
between the two coupled qubit spins. It is defined as^ 1 - 



C12 = max{Ai - A 2 - A3 - A 4 , 0}, 



(52) 



where the quantities Ai > A2 > A3 > A4 are the square 
roots of the eigenvalues of the operator 



R 12 = Ps (a* <g> a y ) p* {a y ® a y ) 



(53) 



n=0 



We find Ai are values in decreasing order of ^/pi 1/944 + 
IpwIj \y/pupH— Ipull, 2/922, and 0. For the system with 
an initial state of |-0) = |00), i.e., the both spins in the 
ground state, we plot the time evolution of the concur- 
rence in Fig. 1. Although there is no initial entanglement 
and no coupling between the two spins, it is interesting to 
notice that the entanglement between the two spins after 
some time is present as shown in Fig. 1(a). This confirms 
that the environment which usually causes the decoher- 
ence of the system can nevertheless entangle qubits that 
are initially prepared in a separable stat o 39 i 42 ' 43 i 44 ' 45 . 
This is mainly due to the fact that the two spins are 
coupled to the same, common environment which then 
in turn generates some effective interaction between the 
two spins even if they were originally decoupled. The 
result, however, depends on the environmental temper- 
ature. Further numerical calculations show that no en- 
tanglement is generated, for example, for T > 8g. As the 
coupling between the two qubit spins is switched on even 
though the value is small, the "collapse" and "revival" of 
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FIG. 1: (color online) Time evolution of concurrence for an 
initial two-qubit state of \tp) — 1 00) for different values of Q; 
(a) Q, — and (b) f2 = 0.03g. Other parameters are fio = 2g, 
90 = g, T = lg. 



the entanglement as a function of time are demonstrated 
in Fig. 1(b). This is in analogy to the collapse and revival 
of atomic population inversion of a single two-level atom 
interacting with a single mode field initially in a coher- 
ent stated, a Fock stated or a squeezed stated in quan- 
tum optics. Here from the Hamiltonian, Eqs. (2), (10) 
and (11), this novel phenomenon of entanglement arises 
from two-coupled qubits interacting with a single mode 
field initially however in a thermal state. The reasons for 
causing the collapse and revival behaviors in these cases 
could be similar, that is the Rabi (or time evolution) 
oscillations associated with different excitations have dif- 
ferent frequencies . Consequently, as the time increases, 
these Rabi (or time evolution) oscillations become uncor- 
rected leading to a collapse behavior. As time is further 
increased, the correlation is restored and the revival oc- 
curs. 

Figure 2 shows the time evolution of the concurrence 
for the system in the initial state of \if>) = -^=(|00) + |11)), 
i.e., an maximally entangled state. At high temperature, 



the state loses its entanglement completely for a short 
period of time, and then it is partially entangled again 
some times later. This is in agreement with the results 
of Ref. GJ where an initially unpolarized infinitely tem- 
perature states of the spin bath is assumed. However, 
our results are temperature dependent. At a very low 
temperature, the concurrence exhibits the behavior of 
oscillation between 1 and 0.35. With increasing temper- 
atures, the concurrence decreases more quickly and os- 
cillates disorderly in the lower value region. At a fixed 
time t, the concurrence decreases with temperatures and 
a critical temperature T c exists, above which the entan- 
glement vanishes. However, this critical temperature T c 
is time-dependent and sensitive to the initial state of the 
system. Figure 3 illustrates the time evolutions of the 
concurrence for different values of the coupling constant 
go- As expected, increasing the value of the coupling 
constant has similar effects as increasing the value of 
the environmental temperature, i.e., the decay rate of 
the concurrence increases. Figure 4 presents the time 
evolution of the concurrence for different inner-bath-spin 
coupling constants g. We see that the concurrence in- 
creases with the increase of g. This confirms that strong 
quantum correlations within the environment suppress 
decoherenc o 24 ' 48 i 49 and thus perhaps also disentangle- 
ment. As shown in the inset, the concurrence is regained 
sometime later. However, it appears disorderly without 
a particular pattern. Similar behaviors arise for Fig. 2, 
Fig. 3, and Fig. 5 in the long time scales, which reflect 
the Non-Markovian dynamics of the system. In Fig. 5, 
we show the effect of the coupling between the two qubit 
spins on the concurrence. It is obvious that the coupling 
benefits the entanglement. Note that the initial state of 
the system in this case is different from that in Fig. 1. 
So, if the system is initially prepared in a maximally en- 
tangled state, the larger the coupling constant £1 is, the 
more slowly the entanglement decays. 

When the two spins are initially prepared in their ex- 
cited state, i.e. \tp) = |11), the result is quiet different 
from that in Fig. 1. If the detuning (1q = 2<?, i.e., at the 
on-resonant case, no entanglement between the two qubit 
spins exists at any temperatures, even with a strong in- 
teraction fl between them. This is consistent with the re- 
sult for two qubit atoms obtained by Kim et al. ^ where 
the coupling between the two qubit atoms is not consid- 
ered and the detuning between the two atoms and the 
single-mode field is zero. If the detuning fiQ 2g, i.e., 
in the off-resonant case, the two qubit spins will entangle 
via the environment again as shown in Fig. 6. So the 
entanglement generation of the two spins in this case is 
very sensitive to the detuning which can be controlled in 
our model by the locally applied external magnetic field. 



IV. CONCLUSION 

We have studied the exact entanglement evolution of 
two coupled qubit spins in a model of a quantum Hciscn- 
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FIG. 2: (color online) Time evolution of concurrence for an 
initial two-qubit state of \tjj) = ( |00) + |H}) for different 
temperatures; T = lOg (solid curve), T = 5g (dashed curve) 
and T = O.lg (dot dashed curve). Other parameters are fio — 
2#, go = g, O = 0. 
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FIG. 3: (color online) Time evolution of concurrence for an 
initial two-qubit state of \ip) — -^=(|00) + |11)) for different 
values of go; go = g (solid curve), go = 0.5g (dashed curve) 
and go = 0.2g (dot dashed curve). Other parameters are 
Ho = 2g,T = 5g, Q, = 0. 



berg XY spin star environment in the thermodynamic 
limit. The dynamics of the reduced density matrix of the 
two coupled spins is analytically obtained in terms of a 
novel operator technique which is mathematically simple 
and physically clear. In our analysis, the transformed 
Hamiltonian becomes effectively Jaynes-Cumming like 
and thus the results are also relevant to cavity quantum 
electrodynamics. 

The time evolutions of the concurrence of the two cou- 
pled spin qubits for different initial conditions are eval- 
uated exactly. The results show that the dynamics of 
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FIG. 4: (color online) Time evolution of concurrence for an 
initial two-qubit state of \ip) = -i=(|00) + |11)) for different 
values of g; g — go (solid curve), g = 1.5po (dashed curve) and 
g = 2go (dot dashed curve). Other parameters are /xo = 2g, 
T — 5go, fi = 0. The inset shows the long time behavior of 
concurrence for g = go- 
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FIG. 5: (color online) Time evolution of concurrence for an 
initial two-qubit state of |i/>) = -j=(|00) + |11)) for different 
values of Q; = (solid curve), Q = 3g (dashed curve), and 
SI = 6g (dot dashed curve). Other parameters are fio = 2g, 
go =g,T= 5g. 



the entanglement strongly depends on the initial state of 
the system, the coupling between the two spin qubits, 
the interaction between the qubit system and the en- 
vironment, the interactions between the constituents of 
the spin environment, the environmental temperatures, 
as well as the detuning controlled by a locally applied ex- 
ternal magnetic field. We have also found that if the two 
coupled spin qubits are initially prepared in the ground 
state, the entanglement between them will exhibit the 
"collapse" and "revival" behavior with time due to the 
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FIG. 6: (color online) Time evolution of concurrence for an 
initial two-qubit state of — |11) for different values of de- 
tuning; /io = 3.5g (solid curve), and fio = Gg (dashed curve). 
Other parameters are go — g, T — lg, £7 = 0. 



interaction between the two spin qubits and the environ- 
ment, in analogy to the collapse and revival of the atomic 
population inversion in quantum optics. 
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